Claims : 

1. A method of processing a digital image, comprising the 
steps of: 

providing digital data indexed to represent positions of 
an image having S spectral bands for simultaneous output on a 
display, said digital data being indicative of an intensity 
value l^Xjy) for each position (x,y) in each i-th spectral 
band ; 

defining a classification of said image based on dynamic 
range of said image in each of said S spectral bands; 

adjusting said intensity value for said each position in 
each i-th spectral band to generate an adjusted intensity 
value for said each position in each i-th spectral band in 
accordance with 

N 

Y,W n (log l ± (x,y) - log [J^y) *F n {x f y) ] ) , i = i,..,s 

n=l 

where S is the number of unique spectral bands included in 
said digital data and, for each n, W n is a weighting factor 
and F n (x,y) is a unique surround function applied to said each 
position (x,y) and N is the total number of unique surround 
functions ; and 

filtering said adjusted intensity value for said each 
position of said image in each of said S spectral bands using 
a filter function based on said classification of said image 
wherein a filtered intensity value R f (x,y) is defined. 
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2 . A method according to claim 1 wherein each said unique 
surround function is a Gaussian function. 



1 3. A method according to claim 2 wherein said Gaussian 

2 function is of the form 

3 satisfying the relationship 

-i 2 

k n f f e Cn dxdy = 1 



B *4 where 

35 i = \Jx 2 +y 2 

and, for each n, k n is a normalization constant and c p is a 
■■*£> unique constant for each of said N unique surround functions. 

h1 4. A method according to claim 1 further comprising the step 

12 of multiplying said filtered intensity value Rj (x,y) by 

[ BI ± (x,y) 
log 

E JiU # y) 

3 to define a color-restored intensity value R ! -(x f y), where B 

4 is a constant. 

1 5. A method according to claim 1 wherein said each position 

2 (X/Y) defines a pixel of said display. 
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6. A method according to claim 1 wherein, for each n, W =1/N. 

7. A method according to claim 1 wherein said step of 
defining comprises the step of using image statistics 
associated with said image in each of said S spectral bands to 
select said filter function. 

8. A method according to claim 7 wherein said image 
statistics include brightness and contrast of said image in 
each of said S spectral bands. 

9. A method according to claim 1 further comprising the steps 
of: 

selecting a maximum intensity value V-(x,y) from the 
group consisting of said intensity value I-(x,y) and said 
filtered intensity value Rj (x,y); and 

displaying an improved image using said maximum intensity 
value V 1 (x , y ) . 

10. A method according to. claim 4 further comprising the 
steps of: 

selecting a maximum intensity value V^x^y) from the 
group consisting of said intensity value I f (x,y) and said 
color-restored intensity value R ! -(x,y); and 

displaying an improved image using said maximum intensity 
value V i (x, y) . 



11. A method of processing a digital image, comprising the 
steps of: 

providing digital data indexed to represent the positions 
of a plurality of pixels of a J-row by K-column display, said 
digital data being indicative of an intensity value l(x,y) for 
each of said plurality of pixels where x is an index of a 
position in the J-th row of said display and y is an index of 
a position in the K-th column of said display wherein a JxK 
image is defined; 

convolving said digital data associated with each of said 
plurality of pixels with a function 

-r 2 
r 2 

e c 

to form a discrete convolution value for each of said 
plurality of pixels, said function satisfying the relationship 



k is a normalization constant and c is a constant; 

converting, for each of said plurality of pixels, said 
discrete convolution value into the logarithm domain; 

converting, for each of said plurality of pixels, said 
intensity value into the logarithm domain; 



kite c dxdy = 1 




where 



r = 




20 subtracting, for each of said plurality of pixels, said 

21 discrete convolution value so-converted into the logarithm 

22 domain from said intensity value so-converted into the 

23 logarithm domain, wherein an adjusted intensity value is 

24 generated for each of said plurality of pixels; and 

2 5 filtering said adjusted intensity value for each of said 

2 6 plurality of pixels with a filter function that is based on 

27 dynamic range of said JxK image wherein a filtered intensity 

28 value R(x,y) is defined. 



12. A method according to claim 11 wherein the value of said 
SI constant c is selected to be in the range of approximately 

^ 0.01 to approximately 0.5 of the larger of J and K, 

f3| 13 . A method according to claim 11 further comprising the 
steps of: 

JSf selecting, for each of said plurality of pixels, a 

f t maximum intensity value V(x,y) from the group consisting of 

5 said intensity value I(x,y) and said filtered intensity value 

6 R(x,y) ; and 

7 displaying an improved image using said maximum intensity 

8 value V(x,y) . 
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1 14 . A method of processing a digital image, comprising the 

2 steps of: 

3 providing digital data indexed to represent the positions 

4 of a plurality of pixels of an J-row by K-column display, said 

5 digital data being indicative of an intensity value I f (x,y) 

6 for each i-th spectral band of S spectral bands for each of 

7 said plurality of pixels where x is an index of a position in 

8 the J-th row of said display and y is an index of a position 

9 in the K-th column of said display wherein a ( JxK) f image is 
10 defined for each of said S spectral bands and a JxK image is 
LI defined across all of said S spectral bands; 

|f defining a classification of said JxK image based on 

|f dynamic range of each said (JxK) - ; 

14 convolving said digital data associated with each of said 

15 plurality of pixels in each i-th spectral band with a function 



for n-2 to N to form N convolution values for each of said 

17 plurality of pixels in each said i-th spectral band, said 

18 function satisfying the relationship 

k n ffe °" dxdy = 1 

19 where 

i = \jx 2 +y 2 

20 and, for each n, k n is a normalization constant and c n is a 



25 



21 unique constant; 

22 converting, for each of said plurality of pixels in each 

23 said i-th spectral band, each of said N convolution values 

24 into the logarithm domain; 

25 converting, for each of said plurality of pixels in each 
2 6 said i-th spectral band, said intensity value into the 

2 7 logarithm domain; 

28 subtracting, for each of said plurality of pixels in each 

29 said i-th spectral band, each of said N convolution values so- 
il converted into the logarithm domain from said intensity value 

:, 

i : J so-converted into the logarithm domain, wherein an adjusted 

|| intensity value is generated for each of said plurality of 

|;i pixels in each said i-th spectral band based on each of said 

3 4 N convolution values; 

forming a weighted sum for each of said plurality of 

5j€ pixels in each said i-th spectral band using said adjusted 

§■=? intensity values; and 

38 filtering said weighted sum for each of said plurality of 

39 pixels in each said i-th spectral band with a filter function 

40 that is based on said classification of said JxK image wherein 

41 a filtered intensity value R-(x,y) is defined. 

1 15. A method according to claim 14 wherein the value for each 

2 said unique constant c n is selected to be in the range of 

3 approximately 0.01 to approximately 0.5 of the larger of J and 

4 K. 
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16. A method according to claim 14 further comprising the 
step of multiplying said filtered intensity value R { (x,y) by 



log 



Bljjx.y) 



1=1 



3 
4 
5 



to define a color-restored intensity value R , 1 -(x / y), where B 
is a constant and S is a whole number greater than or equal to 
2. 



1 17 . A method according to claim 14 wherein said step of 

-.2 defining comprises the step of using image statistics 

'f§ associated with each said ( JxK) - image to select said filter 

: M function. 

18. A method according to claim 17 wherein said image 

rg statistics include brightness and contrast of each said (JxK) 1 

^3 image. 
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1 19. A method according to claim 14 further comprising the 

2 steps of: 

3 selecting a maximum intensity value V-(x,y) from the 

4 group consisting of said intensity value I f (x,y) and said 

5 filtered intensity value R-(x,y) ; and 

6 displaying an improved image using said maximum intensity 

7 value V 1 - (x, y) . 

1 20. A method according to claim 16 further comprising the 

2 steps of: 

y3$ selecting a maximum intensity value V f (x,y) from the 

SI group consisting of said intensity value I f (x,y) and said 

%f> color-restored intensity value R } i (x,y) ; and 

y5 displaying an improved image using said maximum intensity 

H7 value V-(x,y). 
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1 21. A method of processing a digital image, comprising the 

2 steps of: 

3 providing digital data indexed to represent positions of 

4 an image having S spectral bands for simultaneous output on a 

5 display, said digital data being indicative of an intensity 

6 value l^(x,y) for each position (x,y) in each i-th spectral 

7 band ; 

8 defining a classification of said image based on dynamic 

9 range of said image in each of said S spectral bands; 

10 adjusting said intensity value for said each position in 

11 each i-th spectral band to generate an adjusted intensity 
ik value for said each position in each i-th spectral band in 
£B accordance with 

I? N 

f 3 J2 w n < lo 9 ±i(x<y) ~ log[I 2 (x,y) *F n (x,y)] ) , 1 = 1, ,S 

ft where S is a whole number greater than or equal to 2 and 

jt$ defines the total number of spectral bands included in said 

t% digital data and, for each n, W n is a weighting factor and 

17 F n (x,y) is a unique surround function of the form 

-r 2 

e * 

18 satisfying the relationship 



k n J J e Cn dxdy = 1 

19 where 
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r = \]x 2 +y 2 

2 0 and, for each n, k n is a normalization constant and c n is a 

21 unique constant where N is the total number of unique surround 

22 functions; 

2 3 filtering said adjusted intensity value for said each 

24 position in each i-th spectral band with a function based on 

25 said classification of said image wherein a filtered intensity 

26 value R-(x,y) is defined; and 

27 multiplying said filtered intensity value R,- (x,y) by 

f 51,(^7) ' 

U log 

$6 to define a color-restored intensity value R ! -(x,y), where B 

SSS is a constant. 

22. A method according to claim 21 wherein, for each n, 
;;| W n =l/N. 



1 2 3. A method according to claim 21 wherein the value for each 

2 said unique constant c n is selected to be in the range of 

3 approximately 0.01 to approximately 0.5 of the larger of J and 

4 K. 
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1 24. A method according to claim 21 wherein said step of 

2 defining comprises the step of using image statistics 

3 associated with said image in each of said S spectral bands to 

4 select said filter function. 

1 25. A method according to claim 24 wherein said image 

2 statistics include brightness and contrast of said image in 

3 each of said S spectral bands. 

^4. 26 . A method according to claim 2 1 further comprising the 
steps of: 

% selecting a maximum intensity value V f (x,y) from the 

3 group consisting of said intensity value I-(x,y) and said 

V s *5 color-restored intensity value R'-(x,y); and 

f; : ^6 displaying an improved image using said maximum intensity 

/IJ7 value \f-(x,y). 



31 



